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ABSTRACT 

This paper presents a careful and detailed light curve analysis of bright detached 
eclipsing binaries (DEB) in the Small Magellanic Cloud, discovered by OGLE collab- 
oration, on the basis of recently available difference image analysis (DIA) photometry. 
The 19 binaries brighter than 16.4 mag in I band and with the depth of primary 
and secondary eclipse greater than 0.25 mag were investigated. The solutions were 
obtained by a modified version of the Wilson-Devinney program. The quality of DIA 
light curves - a good phase coverage and relatively small scatter - is enough to cal- 
culate realistic estimates for the third light I3 and the argument of periastron lo . It 
was found that solutions of detached, eccentric systems with flat light curve between 
eclipses usually may suffer from indetermination of Z3 in contrast to those of similar 
systems having some proximity effects. 

The physical properties of the stars were estimated on the basis of their pho- 
tometric elements and indices assuming the distance modulus to SMC ~ 18.9 and 
consistency between an empirical mass-luminosity relation and the flux scaling. The 
method was tested on three LMC stars of known absolute dimensions and a good 
agreement was found for m — M ~ 18.5. Such an approach may give fast and accurate 
estimates of absolute dimensions for large and homogeneous samples of eclipsing bina- 
ries in the Magellanic Clouds and other close galaxies. Moreover, this method allows 
also for independent estimation of E{B — V) in the direction to a particular binary. 
The subset of six bright DEB's worth future intensive investigations as likely distance 
indicators to SMC, was chosen. They are SC3 139376, SC4 53898, SC5 129441, SC6 
67221, SC6 215965 and SC9 175336. 
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1 INTRODUCTION 

Accurate determination of absolute dimensions of eclipsing 
binaries based, among others, on detailed light curve anal- 
ysis is still a very important task of modern astronomy. It 
gives an opportunity for testing advanced evolutionary mod- 
els of stars (e.g. Pols et al. 1997) but also allows for very pre- 
cise distance determination. The method, at least in princi- 
ple, is quite straightforward (see Paczyhski 1997 for an out- 
line and Clausen 2000 for a review) . Three distance determi- 
nation for eclipsing binaries in the Large Magellanic Cloud 
were quite recently presented: HV 2274 (Guinan et al. 1998) , 
HV 982 (Fitzpatrick et al. 2001) and EROS 1044 (Ribas et 
al. 2002) giving the distance modulus m — M ~ 18.4. Light 
curve solutions together with estimation of physical param- 
eters and the distance to the Small Magellanic Cloud were 
obtained also for two eclipsing binaries: HV 2226 (Bell et 
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al. 1991) and HV 1620 (Pritchard et al. 1998) giving for the 
distance modulus m — M ~ 18.6. Very recently ten eclipsing 
binaries have been analysed in the central part of the SMC 
giving mean distance modulus m — M = 18.9±0.1 (Harries, 
Hilditch & Howarth 2003). The light curve analysis in their 
paper was based on DoPHOT OGLE photometry (Udalski 
et al. 1998). 

In recent years, extensive CCD photometry have been 
undertaken during the projects of microlensing searching 
toward the Magellanic Clouds and Galactic Bulge. As a 
byproduct, thousands of eclipsing binaries were discovered 
and placed in catalogs (Grison et al. 1995, Alcock et 
al. 1997, Udalski et al. 1998). Preliminary solutions for 
MACHO LMC catalog were found by Alcock et al. (1997). 
They fitted 611 light curves of eclipsing binaries using 
the EBOP code (Etzel 1993). That code uses a tri-axial 
ellipsoid approximation what makes this code very fast. 
But the code cannot properly account for the proximity 
effects present in most of the eclipsing binaries detected by 
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MACHO. For SMC eclipsing binaries OGLE collaboration 
(Udalski et al. 1998) extracted from their catalog a sample 
of 153 detached binaries with well defined, narrow eclipses 
of similar depth and relatively good photometry. It was 
argued that the sample contains the best systems for 
distance determination. Wyithe & Wilson (2001, 2002 - 
hereafter WW1, WW2) presented an excellent analysis of 
the whole OGLE data containing 1459 eclipsing binaries 
using an automated version of the Wilson- Devinney (1971, 
Wilson 1992) code. This program can provide very detailed 
treatment of proximity effects. They used public domain 
DoPHOT photometry available at OGLE homepage 
http :/ / sirius. astrouw. edu.pl/ ogle / ogle2/var- stars /smc/ecl 
in the analysis. They extracted in an objective way two 
subsets of detached binaries with very narrow eclipses 
and complete eclipses as the most likely good distance 
indicators. Subset of semi-detached eclipsing binaries was 
also chosen as distance indicators on the basis of similar 
surface brightness between components and complete 
eclipses. The authors intend, also, to extract such subsets 
for contact binaries. In general it is believed that the 
close-to-ideal system chosen as a distance indicator, should 
have two stars of similar temperature and luminosity and 
the separation of the components should be significantly 
larger than their radii. But well-detached, eccentric systems 
may suffer from the aliasing of photometric solutions. WW2 
argued that semi-detached and contact binaries, and also 
detached eclipsing binaries showing larger proximity effects, 
have their photometric solutions better defined and may 
serve as potential best distance indicators as in the case 
of recently analyzed 15-th magnitude binary EROS 1044 
(Ribas et al. 2002). In this paper I focused my attention 
only on detached eclipsing binaries. 

The release of new DIA photometry by OGLE collabo- 
ration motivated me to reanalyse the light curves of eclips- 
ing binaries in SMC for the purpose of extracting a sample 
of the best distance indicators. This sample was chosen us- 
ing additional important criteria which were not used in 
the selections made by Wyithe & Wilson. The analysis was 
restricted only for bright systems. It is obvious that those 
systems have, in general, the largest probability of obtain- 
ing high-quality light curves and radial velocity curves. Thus 
the searching for distance indicators should start from the 
analysis of the brightest candidates. Moreover using DIA 
photometry we can achieve more realistic determination of 
the light curve parameters than in the case of using the older 
DoPHOT photometry. 



2 THE DATA 
2.1 Photometry 

CCD differential photometric observations of SMC eclips- 
ing binaries were obtained by OGLE between 1997 
and 2000. The data were acquired using 1.3 m War- 
saw Telescope at Las Campanas Observatory. The pho- 
tometric / indices were calculated using the technique 
of difference image analysis (DIA) as was reported in 
Zebruh, Soszyriski & Wozniak (2001a). The new pho- 
tometry was quite recently facilitated on OGLE home- 
page http://sirius.astrouw.edu.pl/ogle/ogle2/dia (Zebruh et 
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Figure 1. The light curves of SMC SCI 25589 eclipsing bi- 
nary based on OGLE DoPHOT photometry (upper panel) and 
OGLE DIA photometry (lower panel). Preliminary solutions were 
marked by solid lines. 



al. 2001b). The OGLE database contains all variables de- 
tected in the LMC and SMC (68 thousand) during the course 
of the OGLE-II project. The DIA is argued as a superior 
method over classical DoPHOT method because of its better 
accuracy, especially in crowded fields (Zebruh et al. 2001a). 
The light curves based on OGLE-DIA photometry are char- 
acterized by considerably smaller scatter of observations and 
higher smoothness of the light curves. Also, in most cases 
the phase coverage is better because the DIA photometry 
includes OGLE observations which have extended almost 
two years after the moment when original catalog of SMC 
eclipsing binaries (Udalski et al. 1998) was published. 

The quality of OGLE-DIA photometry can be demon- 
strated by a comparison of the / light curves obtained with 
both methods. As an example I present a 17-th magnitude 
eccentric eclipsing binary SMC SCI 25589 - Fig. 0- together 
with preliminary light curve solutions. Note that according 
to the analysis of DIA light curve the system has apparently 
total eclipses and much larger separation between compo- 
nents than in the case of the analysis of DoPHOT data. 

The candidates for the analysis were chosen after the 
visual inspection of the light curves of eclipsing binaries in 
SMC released by Udalski et al. (1998). A preliminary set 
of about 25 binaries were selected brighter than / = 16.4 
mag at maximum brightness and with the depth of the pri- 
mary and secondary eclipses greater than 0.25 mag. The 
candidates fall into three separate groups. The first group, 
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Figure 2. The light curves of SMC SC6 67221: from Udalski 
et al. (1998) (upper panel) and from the analysis of OGLE-DIA 
photometry in this paper (lower panel). 



hereafter called A, consist of well-detached, eccentric sys- 
tems with, usually, deep narrow eclipses. The second group, 
hereafter called B, comprise the short-period systems with 
apparently circular orbits and the small separation between 
components. As a result, the systems from group B show 
quite large proximity effects in their light curves. The third 
group C comprises well separated long period systems with 
circular orbits. As we expect long period systems to be, in 
general, eccentric this is most probably a chance selection. 

Afterwards, the photometric data for each selected star 
were extracted, using the coordinates as identifier. Some 
stars were not found in OGLE-DIA database because they 
probably passed undetected through the variability check- 
ing filters (I. Soszyriski 2002, private communication). The 
Table^gives the names and identification of 19 selected sys- 
tems together with a summary of photometric parameters. 

2.2 Initial data preparation 

The observations were folded with trial periods following the 
method given by Kaluzny et al. (1998). The searches were 
usually restricted to a period close to the period given by 
Udalski et al. (1998). The final period was found by tuning 
to the shape of the minima. It was found that in most cases 
the resulting periods agreed very well with those given by 
Udalski et al. (1998). However, for SC6 67221, I have found 
a new period resulting from observations of both eclipses. 
This binary turned out to be a very eccentric system with 
the orbital period of 10.977 days. The period is close to a 
whole number of 11 days and this is the reason for spurious 
secondary eclipse detection in the original OGLE catalog. 
Fig. [5] presents the light curve based on the proper period 
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Figure 3. The light curves of SC5 311566 tuned to the shape of 
primary minimum (upper panel) and to the secondary minimum 
(middle panel). Lower panel: the light curve computed assuming 
a progressive shift of the secondary minimum in respect to the 
primary minimum - see the text for explanation. 



and, for comparison, the light curve from the original OGLE 
catalog. 

Two eccentric binaries SC3 139376 and SC5 311566 
show small changes in the light curve most probably due 
to the apsidal motion. In both cases the tuning to the shape 
of the primary and the secondary minimum gives signifi- 
cantly different periods - see as an example SC5 311566 
(FigEJ. Smooth light curves were obtained assuming that 
the change of the shape of minima is small and the sec- 
ondary minimum was progressively shifting in respect to the 
primary minimum during the period of observations. The 
linear correction was applied to all phases corresponding to 
the secondary eclipse (phases: 0.35-0.55 for SC3 139376 and 
0.4-0.6 for SC5 311566) in the form: 

<t>i = <t>i + a X E, (1) 
where eft and E are the phase and the epoch number of the 
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Table 1. The selected eclipsing binaries 
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SC3 139376 


0:44:08.69 
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6.05213 


629.07158 


14 
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-0.165 


-0.126 


SC5 129441 


0:49:40.56 


-73:00:22.5 


8.05104 


474.20881 


16 
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-0.105 


-0.072 


SC5 311566 


0:51:34.84 


-72:45:45.9 


3.29132 


468.56556 


16 
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-0.120 
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SC6 11143 


0:51:39.69 


-73:18:44.8 


5.72606 


471.95035 


16 


.251 


-0.138 


-0.092 


SC6 67221 


0:52:06.22 
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10.97765 


466.88188 


15 
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-0.101 


-0.095 


SC6 221543 


0:53:40.39 
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3.41680 


468.59872 


16 


.089 


-0.218 


-0.109 


SC6 272453 


0:54:33.22 
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5.70970 


469.30090 


15 


.805 
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-0.057 


SC9 163575 


1:02:46.27 
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1.97109 


629.49776 


15 
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-0.221 


SC10 37223 


1:03:41.42 
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2.58124 


630.63675 


16 
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The B subset 
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16 
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14 


.867 
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SC4 192903 


0:48:22.69 


-72:48:48.6 


102.86800 


761.00300 


16 


.242 


0.771 


0.956 


SC8 201484 


1:00:18.05 
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14 
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Calculated by the period finding program. See H2.21 
Calculated from light curve analysis. See § 3.1 



z-th observation. Denoting the period resulting from the tun- 
ing to the primary eclipse by P p and to the secondary one by 
P s , the factor a was calculated by formula: a = (P p — P 3 )/ P p 



and was 7.5 ■ 1(F 5 and 7.0 
311566, respectively. 



1CT J for SC3 139376 and SC5 



The old reference epochs T a given by Udalski et 
al. (1998) were used to fold up the observations. New ref- 
erence epochs were determined later during the light curve 
analysis - see § 3.1. For each star I calculated out-of-eclipse I 
magnitude at the maximum of the brightness. The DIA data 
does not include the B and V photometry. For the purpose of 
obtaining the most accurate B — V and V — I colours I recal- 
culated the out-of-eclipse B, V and I magnitudes using the 
old DoPHOT photometry. Note, that the magnitudes which 
I have found refer to the mean magnitudes at the maximum 
brightness rather than to the maximal magnitudes given in 
the original OGLE catalog. Resulting out-of-eclipse / mag- 
nitudes were then used to normalize the light curves. 

In two cases (SC6 11143 and SC9 163575) the light 
curve shows considerable amount of scatter, much larger 
than the average error of the observations. The light curves 
of both systems are presented on Fig. 2] (upper panel). I have 
removed observational points when the eclipses occurred and 
rerun the period finding program to find some periodicity in 
the scatter. It turns out that in the case of SC6 11143 the 
out-of-eclipse observations assemble with the period of 2.618 
days and form a light curve with two well-defined, narrow 
minima of different depth (middle panel of Fig.0J. It is likely 
that the binary SC6 11143 blends with another short period 
and slightly eccentric eclipsing binary. The similar analy- 



sis done for the system SC9 163575 shows that the out-of- 
eclipse observations assemble with the period of 1.985 days 
and form two minima separated by a half of the period. It 
seems that the system SC9 163575 blends also with a short 
period eclipsing binary. Their periods are very close to each 
other and in the original light curve the points, correspond- 
ing to the minima of the blend, form two clumps situated 
around the minima of SC9 163575. The light curves of both 
systems were cleaned by removing from the light curve all 
points corresponding to the moments when the eclipses oc- 
curred in the blended binary. The resulting light curves are 
presented in the lower panel of Fig. 2] Also the system SC5 
38089 in spite of its relatively high brightness (I ma x = 15.26) 
has quite large scatter of observations. The analysis done for 
this system shows that the changes of the brightness are ran- 
domly distributed and thus, it is likely that the reason of the 
scatter is some kind of an intrinsic variability. 



3 LIGHT CURVE ANALYSIS 
3.1 Method 

All fits to the light curves were obtained using a modified 
version of the Wilson-Devinney(1971, Wilson 1979, 1992, 
hereafter WD) program with a model atmosphere proce- 
dure based on the ATLAS9 code and developed by Milone, 
Stagg & Kurucz (1992). An automated iterations scheme of 
WD program was employed. For A subset a detached con- 
figuration (mode 2) was chosen for all solutions and simple 
reflection treatment was used (MREF=1, NREF=1). In the 
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Figure 4. Upper panel: the original DIA light curves of SC6 11143 and SC9 163575. Middle panel: the out-of-eclipse observations of 
both stars folded up with a period and an epoch given in each picture. Both light curves suggest that a blend is an eclipsing binary. 
Lower panel: the final light curves after removing the variability of the blend. 



case of B subset both detached and semi-detached (mode 5) 
configurations were used to obtain solutions. For B and C 
subsets the detailed reflection treatment with a double re- 
flection was taken into account (MREF=2, NREF=2). For 
all stars logg = 4.0 was fixed except for stars from C subset 
for which logp = 2.0 was assumed. 

I have found the initial guess via trial-and-error fitting 
of the light curves. The convergence is defined to have been 
achieved if the parameters corrections given by WD pro- 
gram are smaller than 0.75 times the standard errors (1.0 
in the case of eccentricity and argument of periastron) on 
two consecutive iterations. The adjusted parameters were 
the surface potentials (fii, f^), the effective temperature of 
the secondary (T2), the luminosity of the primary (I/i), or- 
bital inclination (i) and the zero epoch offset (A0). For A 
subset I have adjusted two more parameters: eccentricity (e) 
and the argument of periastron (lo ). 

If for an eccentric system the convergence was not 



achieved I suppressed the adjustment of e, lo and A<f> pro- 
vided the correction were smaller than the standard error for 
these parameters on three consecutive iterations. It turned 
out that such procedure allowed for better convergence and 
the solution was achieved in all cases. Afterwards I have run 
additional iterations, allowing all parameters to vary, to find 
error estimates of the final solution. Similar procedure was 
employed for the systems from B and C subsets, just fixing 
the circular orbit and the suppressing included only the zero 
phase offset. Among the systems in B subset the third light 
(h) adjustment was needed to get the convergence in most 
cases. For all systems correction AT to the epoch given by 
Udalski et al. (1998) was applied in the form: AT = A(j>- P, 
where A<j> is the zero epoch offset of the final solution and 
P is the period of the system. 

WW1 introduced a quantity F e which is greater than 
unity for systems with complete eclipses: 
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r g +r s -cosi 

Fe = to s ' (2) 

where r g is the fractional the radius of larger star and r s is 
the fractional radius of the smaller star, both expressed in 
units of semi-major axis A. F e is valid, in the form presented 
above, for circular orbits or perpendicular orbits (z=90°) 
and also for eccentric orbits provided the fractional radii of 
components are expressed in units of the distance between 
stars at the moment of the mid-eclipse. In general case of 
eccentric orbits, if one eclipse occurs in the neighbourhood 
of periastron and is complete, the other may be partial and 
in the extreme cases, there may be an eclipse near perias- 
tron and none near apastron. The instantaneous distance 
between components d, in units of semi-major axis A, is 
given by the equation: 



1 + e cos v 

where v is the true anomaly. The primary and secondary 
minima occur when the true anomalies v\ and V2 are (ne- 
glecting higher powers of e, see Kopal 1959, Chapter 6, equa- 
tion (9-22)): 

v\ = 90° — to — e cos lj cot 2 i (1 — e sin lo csc 2 i), (4) 
V2 = 270° — lu + e cosoj cot 2 i (1 — e smw csc 2 i), (5) 

respectively. Substituting v\ and v-z into Eq. EJthe distances 
were calculated at the moments of both mid-eclipses. Then 
appropriate quantities and F 2 referring to the primary 
and the secondary eclipse, respectively, were computed. 

Semidetached configuration (mode 5) were investigated 
for all stars from the B subset in order to check which sys- 
tems are definitively separated. When the semidetached con- 
figuration was investigated the mass ratio (q) of the compo- 
nents was adjusted and secondary surface potential (f2a) was 
set by q. However, it turned out that no solution could be 
found in this mode for any star from B subset in spite of 
many tests for different sets of adjustable parameters. Ac- 
cordingly, in subsequent part of the paper I will refer only 
to detached solutions found in mode 2. 



3.2 Parameter choice 

In the analysis presented here I have included the tem- 
perature of binary's components resulting from an intrinsic 
colour (B—V)o of the system. All (B— V) colours from Ta- 
ble Q were dereddened using the average reddening toward 
SMC E(B -V) = 0.087 (Massey et al. 1995, Udalski 2000). 
I set the temperature of the primary (Ti) respectively to the 
(B — V) colour of the system using Flower's (1996) empiri- 
cal calibration. Almost all the selected binaries have minima 
of similar depth what indicates that the difference between 
the temperatures of the components is small. Thus such a 
procedure is justifiable. For those stars for which (B ~ V)o 
colour suggested a very high temperature (~ 4 10 4 K for 
SC5 38089, SC6 215965, SC9 163575, SC9 175336) I set up 
Ti = 35 000 K. 

I used a logarithmic law for the limb darkening (e.g. van 
Hamme 1993). The I limb darkenings 2:1,2, 2/1,2 and bolo- 
metric limb darkenings Xi°%, y\°l were calculated from van 
Hamme's (1993) tables of limb darkenings for appropriate 
log g and T. During the iterations the secondary's limb dark- 



ening coefficients were calculated according to its temper- 
ature T2. The bolometric albedo and gravity brightening 
exponent were set to the canonical value of 1.0 for stars 
with radiative envelopes (von Zeipel 1924). However, for two 
red long-period systems SC4 192903 and SC8 201484 values 
typical for convective envelopes were assumed - bolomet- 
ric albedo was set to 0.5 and gravity brightening exponent 
to 0.32. The synchronous rotation of both components was 
assumed (F\ = F2 = 1). No level depending weighting of 
observations was applied - parameter NOISE was set to 0. 

The mass ratio was assumed to be q — 7712/mi = 1 in 
all calculations. WW2 pointed out that in the case of de- 
tached eclipsing binary the light curve is almost insensitive 
to the mass ratio parameter. Only if one of the components 
is highly evolved (when its radius R increase to about a half 
of the total separation A) or if the system has a very low 
mass ratio, the light curve becomes sensitive for this param- 
eter. The preliminary tests show that no component from 
the set of binaries has fractional radius n = Ri/a larger 
than 0.4. However, for two systems from B subset it was 
necessary to adjust the mass ratio in order to get a good 
solution (see §3.4). 

3.3 Solutions for A subset 

For each system from A group the automated fitting scheme 
described in §3.1 was employed. Only in three cases: SC6 
67221, SC6 221543 and SC6 272453 was the convergence 
achieved by adjusting simultaneously all free parameters. 
The Figure |S] shows the best-fitting model light curves, to- 
gether with O — C residuals, whose photometric parameters 
are given in Table [5] All arguments of the periastron refer 
to the epoch 1999.0. 

For all systems from the A subset reasonable solutions 
were found assuming I3 = 0. However, in the case of SC6 
11143 and SC9 163575 we expect that some amount of the 
third light is present in the system because of the blending 
with another eclipsing variable. In order to estimate third 
light component an additional grid of solutions for different 
values of I3 was calculated in the range of 0.05-0.40 with a 
step of 0.05. For SC6 11143 the best solution was found at 
l 3 = 0.30, while for SC9 163575 - at l 3 = 0.10. It is worth 
noting that the solutions found with I3 = for these systems 
have almost no systematic trends in O-C's residuals, thus the 
intrinsic variability of the blend is the only sign of the third 
light component. 

The solutions for SC6 67221 converged very rapidly and 
had the smallest residuals. Only slight systematic residuals 
are visible near phase 0.04 (just after the egress from pri- 
mary minimum), most probably due to the strong mutual 
reflection effect. Similar but even stronger systematic effect 
is visible near first quadrature ((f) ~ 0.11) in SC5 129441. Be- 
cause of the presence of the proximity effects both systems 
were additionally investigated in order to estimate possibly 
the third light contribution by adjusting I3 as a free pa- 
rameter. In each case synthetic light curves produce a tiny 
improvement only for unphysical values Z3 <0, but such solu- 
tions were invariably unstable. The conclusion is that these 
binaries show virtually insignificant third light contribution 
what is consistent with I3 = 0. For a brief discussion of both 
systems see § 5.1.2. 

Small residuals and relatively quick convergence charac- 
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Table 2. Photometric parameters 
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light ?3 


ratio I2/I1 


pi p2 
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The A 


subset 












SC3 139376* 


0.177±0.004 


0.194±0.004 


1.007±0.006 


83.3±0.2 


0.064±0.003 


133.7±5.7 


0.0 


1 


.191 


0.734 0.704 


SC5 129441* t 


0.158±0.001 


0.162±0.002 


0.999±0.001 


85.2±1.0 


0.375±0.001 


172.0±0.4 


0.0* 


1 


.079 


0.797 0.773 


SC5 311566 


0.157±0.004 


0.147±0.002 


0.962±0.005 


85.7±0.6 


0.251±0.004 


070.0±0.1 


0.0 





.843 


0.841 0.721 


SC6 11143 1 " 


0.133±0.004 


0.112±0.002 


0.972±0.020 


88.8±0.1 


0.263±0.001 


234.3±0.6 


0.30±0.05 


0. 


.704 


0.983 1.022 


SC6 67221* 


0.131±0.002 


0.119±0.002 


1.154±0.011 


83.9±0.3 


0.401±0.002 


043.5±0.9 


0.0* 


1 


.058 


0.757 0.533 


SC6 221543 


0.202±0.006 


0.182±0.005 


0.960±0.007 


84.2±0.7 


0.099±0.002 


146.1±2.7 


0.0 


0. 


.745 


0.794 0.764 


SC6 272453 


0.149±0.004 


0.180±0.004 


0.961±0.006 


85.2±0.6 


0.035±0.004 


236.0±4.5 


0.0 


1 


.392 


0.815 0.831 


SC9 163575 


0.234±0.004 


0.225±0.004 


0.891±0.007 


80.6±0.4 


0.043±0.002 


010.7±7.4 


0.10±0.05 





.762 


0.661 0.655 


SC10 37223 


0.202±0.004 


0.195±0.004 


0.951±0.006 


83.8±0.3 


0.087±0.003 


142.3±2.4 


0.0 


0. 


.874 


0.757 0.728 
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SC4 53898* 


0.189±0.001 


0.356±0.002 


0.910±0.006 


85.7±1.2 


0.0 


090.0 


0.187±0.012 


3 


.035 


1.210 


SC4 103706+ 


0.350±0.002 


0.284±0.006 


0.964±0.004 


86.8±2.0 


0.0 


090.0 


0.382±0.014 


0. 


.612 


1.003 


SC4 163552* 


0.338±0.006 


0.326±0.008 


1.004±0.008 


85.4±0.9 


0.0 


090.0 


0.263±0.007 


0. 


.951 


0.888 


SC5 38089 


0.246±0.009 


0.255±0.007 


0.988±0.003 


76.8±1.9 


0.0 


090.0 


0.0 


1 


.053 


0.545 


SC6 215965* t 


0.277±0.003 


0.335±0.002 


0.922±0.001 


75.6±0.3 


0.0 


090.0 


0.0* 


1 


.356 


0.656 


SC6 319966+ 


0.339±0.006 


0.324±0.008 


1.007±0.007 


88.3±3.0 


0.0 


090.0 


0.462±0.005 





.861 


0.973 


SC8 104222 


0.311±0.011 


0.331±0.009 


0.993±0.008 


86.8±1.4 


0.0 


090.0 


0.087±0.013 


1 


.129 


0.894 


SC9 175336* 


0.194±0.002 


0.358±0.002 


0.965±0.005 


88.7±0.8 


0.0 


090.0 


0.174±0.010 


3 


.228 


1.330 










The C 


subset 












SC4 192903 


0.111±0.004 


0.145±0.003 


0.927±0.004 


87.2±0.7 


0.0 


090.0 


0.0 


1. 


.165 


0.714 


SC8 201484 


0.171±0.003 


0.224±0.002 


0.888±0.002 


78.6±0.2 


0.0 


090.0 


0.0 


1. 


.107 


0.569 



The stars chosen in this paper as likely distance indicators are marked by an asterisk. The stars indicated by Wyithe & Wilson (2001, 
2002) are marked by a dagger. A double dagger means that the binary has most probably no third light contribution. 



terize also SC3 139376 and SC9 163575, whereas, the largest 
residuals are visible in SC5 221543 mostly due to a big in- 
trinsic scatter in observations. 

Having calculated e and 10 we could estimate the rate of 
apsidal motion for two binaries from the sample. The change 
of displacement of minima A$ of an eccentric binary with 
central eclipses (i=90°) is given by (see e.g. Kopal 1959, 
Chapter 6, equations (9-8), (9-9)): 



A$ = 



e sin lo lo 



VT 



(6) 



However for our purpose we can still use this formula also 
for inclination close to 90° provided the timescale of a dis- 
placement change is small with respect to an apsidal period 
U. The change of displacement per one orbital revolution 
(P) is simply a (Eq.0 so the rate of advance of periastron 
is approximated by: 



365 a\/l — e 2 



(7) 



P e sin lo 

Substituting numerical values we get 5?8/yr and l?8/yr for 
SC3 139376 and SC5 311566, respectively. 

3.4 Solutions for B and C subsets 

The photometric parameters for the stars from B and C 
subsets are presented in Table [5] Appropriate light curve 
solutions are shown in Figs. and Q In contrast to the 
A subset, almost all solutions from B subset (beside SC6 
215956) were found only by adjusting the third light Z3 in 
the system. There is a simple explanation of this effect - 
see a discussion in section 16. II SC6 215956 was also inves- 
tigated for possible third light contribution but values of 



I3 suggested by WD program were invariably negative. SC6 
215965 has very small residuals and the solution for this sys- 
tem converged rapidly. Also it seemed to be most stable of 
all in B subset. 

The solution for SC5 38089 was obtained using the spec- 
troscopic mass ratio q = 1.12 and more appropriate tem- 
peratures (Ti,2 ~ 29000 K) determined for this binary by 
Harries et al. (2003). Third light contribution in this system 
is spurious although some improvements to the fits can be 
obtained for £3 > 0.2. 

In spite of the small residuals in SC9 175336 its so- 
lution was extremely unstable and exhibited large system- 
atic O — C deviations, especially near external contacts. The 
tests showed that the reason was the sensitivity of the light 
curve to the mass ratio q as the larger and cooler compo- 
nent seemed to be noticeably more evolved than the hot- 
ter one. The same problems arose for SC4 53898 which 
seems to be similar in many respects to the former, al- 
though its light curve is less sensitive for q. Both systems 
show complete eclipses which allow for more accurate de- 
termination of photometric parameters including the pho- 
tometric mass ratio. Thus I allowed for adjusting of q for 
both systems in detached configuration. The convergence 
was achieved and a stable solution obtained: in the case of 
SC9 175336 for q = 1.15 ±0.05 and in the case of SC4 53898 
for q = 1.36 ±0.12. 

The solutions for two systems from C subset converged 
very rapidly. No third light I3 was adjusted in both cases 
because of the lack of any systematic residuals in the solu- 
tions (Fig.0. Moreover the solution for SC8 201484 has the 
smallest residuals from all of the subsets. 
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Figure 5. The light curves solutions and their residuals for the systems from the A subset. 
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Figure 6. The light curves solutions and their residuals for the systems from the B subset. 
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Figure 7. The light curves solutions and their residuals for the systems from the C subset. 



4 PHYSICAL PARAMETERS 

We cannot directly calculate absolute dimensions and 
masses without knowing spectroscopic orbit. However it is 
possible to estimate stellar parameters using the inversion of 
the method of parallaxes of eclipsing binaries. The method 
of parallaxes was elaborated by Dworak (1974) who based it 
on Gaposchkin (1938) algorithm. We assume that we know 
distance to SMC and we use it to calculate absolute dimen- 
sions by scaling the system to obtain the observed flux at 
Earth. The calculations are quite straightforward. For early 
B or O type stars the Rayleigh- Jeans approximation may be 
used in visible and near-infrared region of the spectrum so 
total V or I luminosity of each star is a linear function of 
temperature Li t 2 ~ Ti,2. Using this approximation we can 
disentangle the temperatures Ti and T2 of both components: 



Ti =T 



l + i&CZVTx)' 



(8) 



where L21 



= L\IL\ = fo/h is the luminosity ratio of 
the components known from the light curve solution (Ta- 
ble |5J and To is the "mean" temperature corresponding 
to the total (B — V)o colour of the system. The temper- 
ature T2 is then directly calculated from the temperature 
ratio. Initially (B — V)o was found for each binary assum- 
ing total, foreground and internal, mean reddening to SMC 
E{B — V) = 0.087 (Massey et al. 1995). The foreground red- 
dening of E(B — V) = 0.037 was assumed (Schlegel, Finke- 
beiner & Davis 1998) with the ratio of the selective to the 
total extinction R — 3.1, while for the internal reddening 
I assumed R — 2.7 (Bouchet et al. 1985). The resulting 
temperatures Ti,2 were used to calculate approximate in- 
trinsic colours (B — V)o and bolometric corrections BC\,2 
of both components via Flower's (1996) calibration for main 
sequence stars. The "standard" mean value m — M — 18.9 
of the distance modulus (DM) for SMC was assumed. The 
semi-major axis A can be then computed by the following 
formula: 



logyl 



0.2 DM - 
0.2 (Mq 



logn - 2 log Ti + 7.524 + 
-Bd- Vi+RE(B-V)), 



(9) 



where Mq is the V absolute magnitude of the Sun (assumed 
+4.75), Vi is the out-of-eclipse V magnitude of the primary 



(calculated assuming h\ jh\ « h/h and that I3 contribu- 
tion in V band is comparable to that in I band). The total 
mass of the system results simply from application of Ke- 
pler's third law. The absolute scale of the system gives im- 
mediately the bolometric luminosity L of each component. 

Now we proceed to the second step of the method. We 
can represent the luminosity of the star as a function of 
its mass, which is well-known as Eddington's relation. This 
relation works only for main sequence stars, but it should 
be kept in mind that for an early spectral type star (massive 
and luminous at the same time) its luminosity is not a strong 
function of the age. The dependence is rather weak and such 
star tends to evolve from ZAMS along almost horizontal 
line on HR diagram. Moreover, no binary from the sample 
(beside two red systems from C subset) seems to contain 
considerably evolved components. 

Fig- OH presents the mass-luminosity relation. The data 
for massive (M > 2.2Mq) Galactic stars were taken from 
compilation of detached, double lined eclipsing binaries done 
by Pols et al. (1997 and references therein) and for three 
LMC eclipsing systems from Ribas et al. (2000), Fitzpatrick 
et al. (2001) and Ribas et al. (2002). The inclusion of 
LMC stars gives an opportunity to check how this rela- 
tion depends on the metallicity. Indeed the Galactic stars 
seems to have a little smaller luminosity then LMC stars 
(having [Fe/H]~ —0.35) at a given mass. For a Galactic 
star the mass-luminosity M — L relation may be approxi- 
mated by a simple linear function in wide range of masses 
(2.2M Q < M < 20M e ): 



logL = 3.664(47) log M + 0.227(38), 



(10) 



i.e. L ~ M 3 ' 66 . The sample of LMC stars do not cover such 
wide range of masses and I assumed that the slope of the 
M — L relation is the same like for the Galactic stars. Then 
the data can be fitted by: 



logL = 3.664(47) logM + 0.380(27) . 



(11) 



Above equation (labeled 2 in Fig. |HJ gives the approximate 
M — L relation for LMC stars. This relation can be com- 
pared with a mass determination for two detached eclipsing 
systems in SMC done by Harries et al. (2003) - see Fig. |S] 
Components of SC6 215965 lie very close to the predicted M- 
L relation, but there are some problems with components of 
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Figure 8. Mass-luminosity relation for Galactic main sequence 
stars (asterisks) given by equation 1 and LMC stars (filled 
squares) — equation 2. Both components of two SMC detached 
systems: SC6 215965 (open squares) and SC5 38089 (open cir- 
cles) from paper Harries et al. (2003) are also marked. Note that 
the components of SC6 215965 are placed almost on the predicted 
M-L relation, while the components of SC5 38089 do not seem to 
concur with this mass-luminosity relationship. 



SC5 38089: especcially a more massive component of 38089 
seems to be a less luminous one. As the physical parameters 
of 215965 are more precisely determined than in the case of 
38089 I decided to adopt the mass-luminosity relation given 
by Eq. llll for SMC stars. Relying on this assumption we can 
check the consistency between parameters found from the 
flux scaling (Eq. |5J and from M — L relation (Eq. 1111 . 

Such treatment allows for independent determination of 
the unknown reddening in the direction to the particular bi- 
nary and thus a reasonable (B — V)o intrinsic colour can be 
obtained. Although the mean foreground reddening toward 
SMC is quite low some large random internal reddenings 
(up to 0.2) were reported for individual B type stars (Larsen, 
Clausen & Storm 2000). It is obvious (e.g. Flower's 1996 cal- 
ibration) that for early type stars the temperature is a very 
strong function of the (B — V)o colour and a small change 
of this colour causes a large change of the temperature. As 
the temperature enters Eq. [5] at a relative high power even 
a small error in determination of reddening E(B — V) pro- 
duces large errors in luminosities and masses. Let us return 
to the problem of the determination of physical parameters. 
We can compute the expected mass ratio q as follows: 

logg= (logL 2 -logLi)/3.664, (12) 

where Li,2 are the absolute luminosities computed by using 
Eq. Then the masses of individual components Mi,2 are 
calculated from Eq. 1111 If there is a consistency between the 
flux scaling and M — L relation a quantity (3: 

(3 = 0.5 (log Li + log L 2 - 3.664 (log Mi + log M 2 )) (13) 

should be equal or close to the free term in Eq. Illl i.e. 0.38. 
I have allowed j3 parameter to vary in the range of 0.36-0.40. 
If f3 were in allowed range the consistency was achieved, the 
E(B — V) = 0.087 was approved as the final reddening and 
calculations stopped. If not, the iterative determination of 
the reddening were done. For j3 greater than 0.40 the correc- 
tion of —0.001 were added to E(B — V), while for /3 smaller 
than 0.36 the correction of +0.001 was applied and all calcu- 



lations were repeated for this new value of reddening, until 
the (3 fell in the allowed range. 

The above method was tested on three LMC eclipsing 
binaries which have their physical parameters known with 
high accuracy. As a free parameter I assumed distance mod- 
ulus to the LMC. The best "fit" to the original data were 
obtained for m — M = 18.5 what is a value still accepted 
as close to the true one (e.g. Gibson 1999) and marginally 
consistent with DM obtained using eclipsing binaries (Ribas 
et al. 2002). The comparison between original data and pa- 
rameters given in Table [3] (lowest panel) show that the "fit" 
is excellent - the method allows for consistent and accu- 
rate estimation of masses, sizes, temperatures, luminosities 
and reddenings of all three systems. It may serve as a proof 
that such an approach of estimation of physical parameters 
works well and may give fast and proper information about 
the absolute dimensions for large samples of eclipsing bina- 
ries - see § 16. 21 for a discussion. It is worth to note that the 
initial value of reddening is irrelevant to the method and the 
same results we obtain adopting another initial E(B — V). 

Table |3] gives the final estimation of the physical pa- 
rameters of each component of binaries in the sample. The 
spectral type and luminosity class of each star was estimated 
from T and log g. SC4 53898 and SC9 175336 were treated in 
a slightly different way. The mass ratio for both systems was 
not calculated from Eq. 1121 but was assumed to be equal to 
the photometric mass ratio determined from the light curve 
solution (see § 3.4). The further calculations were done in a 
similar way like for the rest of the sample. For two binaries 
from C subset which host evolved and cool components the 
reddening was assumed to be equal to the initial value and 
just the simple flux scaling was applied. 

The temperatures Ti from Table |3] in a few cases (SC4 
163552, SC5 311566, SC9 175336) significantly differ from 
initial 21 for which light curve solutions were found. Ad- 
ditional tests show, however, that photometric parameters 
recalculated for the temperatures resulted from the flux scal- 
ing were consistent with those from Table [5] within quoted 
errors. 



5 DISTANCE INDICATORS 

Below the discussion of individual eclipsing binaries as dis- 
tance indicators to SMC is presented. Six candidates were 
chosen according to the following criteria: the presence of 
complete eclipses, small residuals, small and well determined 
third light contribution and the spectral type. Two binaries 
from this subset were also advocated to be possible distance 
indicators by WW1. Another four systems which were se- 
lected by WW1 and WW2 were also briefly discussed in a 
separate subsection, where reasons against their selection 
were given and additionaly SC5 38089 used recently for dis- 
tance determination to SMC is also discussed. 

5.1 Most likely candidates 

5.1.1 SC4 53898 and SC9 175336 

These two binaries may serve as the best candidates be- 
cause of the presence of complete eclipses, high brightness 
(especially SC9 175336), modest proximity effects allowing 
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Table 3. Physical parameters 







Mass 


Mass 


Mean radius 


Temperature 


Luminosity 


Object 


Spectrum 


Mi M 2 

M Q M 


ratio q 


7% log* log 92 


t 2 n 


logLi logL 2 E(B-V) 



SC3 139376* 
SC5 129441* 
SC5 311566 
SC6 11143 1 " 
SC6 67221* 
SC6 221543 
SC6 272453 
SC9 163575 
SC10 37223 

SC4 53898* 
SC4 103706+ 
SC4 163552+ 
SC5 38089 
SC6 215965* t 
SC6 319966+ 
SC8 104222 
SC9 175336* 

SC4 192903 
SC8 201484 

EROS 1044 
HV 982 
HV 2274 



B0V+B0V 
B3III+B3III 
09V+09V 
B2V+B2V 
B2IV+B1IV 
B2V+B2V 
B2V+B2IV 

05+07 
B1V+B1V 

B1V+B2IV 
B0V+B0V 
B1V+B1V 
B0V+B0V 
B0V+B1IV 
B2V+B2V 
B1V+B1V 
B1V+B1IV 



F9IVH 
G0III- 



B2V+B2IV 
B1V+B1V 
B2IV+B2IV 



15.4±2.0 
7.1±0.6 

11.9±1.6 
7.5±0.8 

lO.Oil.l 
8.2±1.0 
7.9±0.9 

16.9±2.4 
8.4±1.0 



16.3±2.1 
7.2±0.6 

11. Oil. 5 
6.6±0.7 

11.1±1.2 
7.4±1.0 
8.4±1.0 

14.6±2.0 
7.8±0.9 



6.5±0.6 8.9±0.9 

10.7±1.3 9.1±1.1 

9.3±1.0 9.2±1.0 

11.7±2.0 H.8±1.9 

14.6±1.1 14.8±1.1 

6.2±0.6 6.1±0.6 

8.4±1.0 8.7±1.1 

lO.Oil.O 11.4±1.1 



G2IV 
G7III 



1.5±0.3 
1.9±0.3 



1.6±0.3 
1.9±0.3 



7.1±0.8 8.7±0.9 
10.4±1.1 10.6±1.1 
12.5±1.3 12. Oil. 2 



1.06±0.11 
1.01±0.07 
0.92±0.10 
0.88±0.08 
1.11±0.09 
0.90±0.09 
1.06±0.09 
0.86±0.10 
0.93±0.08 

1.36±0.12* 

0.86±0.08 

0.98±0.09 

1.01±0.12 

1.02±0.05 

0.98±0.07 

1.03±0.10 

1.15±0.05* 

1.1±0.2 
1.0±0.2 

1.22±0.10 
1.02±0.08 
0.96±0.07 



The A subset 
7.8±0.9 8.6±1.0 
6.5±0.4 6.7±0.4 
4.2±0.5 3.9±0.5 
4.3±0.5 3.6±0.4 
7.5±0.6 6.8±0.6 
4.8±0.5 4.3±0.5 
5.1±0.6 6.1±0.6 
4.9±0.6 4.7±0.6 
4.0±0.4 3.9±0.4 

The B subset 
2.9±0.2 5.4±0.4 
4.9±0.4 4.0±0.4 
5.0±0.5 4.9±0.5 
5.3±0.9 5.5±0.8 
9.0±0.710.9±0.8 
3.4±0.3 3.3±0.3 
4.6±0.6 4.9±0.5 
4.7±0.4 8.7±0.6 

The C subset 
15±2 19±2 
36±2 48±3 
LMC stars 
4.2±0.4 6.7±0.6 
7.0±0.5 7.6±0.5 
10.0±0.7 9.2±0.7 



3.84 


3.79 


31500±1100 


31700±1100 


4.731 


4.823 


0.128 


3.67 


3.65 


16900±400 


16900±400 


3.486 


3.506 


0.094 


4.28 


4.30 


34200±1300 


32900±1200 


4.324 


4.199 


0.180 


4.04 


4.14 


21600±600 


21000±600 


3.565 


3.367 


0.103 


3.69 


3.82 


21700±600 


25100±700 


4.052 


4.218 


0.153 


3.99 


4.03 


22400±600 


21500±600 


3.715 


3.554 


0.027 


3.93 


3.79 


21400±600 


20500±500 


3.680 


3.775 


0.065 


4.29 


4.26 


43300±1800 


38600±1500 


4.875 


4.640 


0.073 


4.15 


4.15 


25100±700 


23900±650 


3.766 


3.648 


0.064 


4.34 


3.93 


24400±700 


22200±600 


3.415 


3.801 


0.122 


4.09 


4.21 


28700±1000 


27700±800 


4.162 


3.917 


0.085 


4.01 


4.03 


24700±700 


24800±700 


3.922 


3.897 


0.216 


4.10 


4.07 


29700±1000 


29400±1000 


4.291 


4.302 


0.047 


3.70 


3.54 


28100±900 


25900±800 


4.653 


4.677 


0.039 


4.16 


4.19 


20700±500 


20800±500 


3.286 


3.259 


0.085 


4.05 


4.01 


23900±650 


23800±650 


3.784 


3.826 


0.113 


4.09 


3.62 


25200±700 


24300±700 


3.901 


4.371 


0.039 


2.26 


2.06 


5900±100 


5400±100 


2.372 


2.472 


0.087 


1.59 


1.36 


5700±100 


5000±100 


3.088 


3.116 


0.087 


4.04 


3.73 


21400±600 


20500±500 


3.524 


3.844 


0.065 


3.77 


3.70 


23300±650 


22800±650 


4.111 


4.143 


0.081 


3.54 


3.60 


23100±650 


23200±650 


4.407 


4.341 


0.121 



Note: all parameters were estimated assuming m — M = 18.9 for SMC stars and m — M = 18.5 for LMC stars. The spectrum was 
estimated from T and log g and, in the case of LMC stars, adopted from the literature. The stars chosen in the present paper as likely 
distance indicators are marked by an asterisk. The stars indicated by Wyithe & Wilson (2001, 2002) are marked by a dagger. A double 
dagger indicates a photometric mass ratio derived from the light curve analysis. 



for reasonable determination of the third light contribution 
to the system and at last, small residuals. In both cases the 
primary is slightly cooler and considerably larger than the 
secondary component. As a result there is a quite high lu- 
minosity ratio between components h/h ~ 3 (see Table |5J 
which may by regarded as the only disadvantage of using 
these systems. The expected velocity semiamplitudes should 
be Ki ~ 250 km/s and K 2 ~ 190 km/s for SC4 53898, while 
Ki ~ 220 km/s and K 2 ~ 190 km/s for SC9 175336. The 
latter lies near the large HII emission region and OB asso- 
ciation - NGC 371. Both stars were not included by WW1 
or WW2 as likely distance indicators. 

5.1.2 SC5 129441 and SC6 67221 

Very high eccentricity and clear proximity effects visible 
near the periastron passage are features of both systems. 
SC5 129441 consists of two somewhat evolved B3 stars of 
equal temperatures. In the second system - SC6 67221 - dur- 
ing the primary eclipse (phase 1.0) the cooler and probably 
larger component is eclipsed and the strong mutual reflec- 
tion causes the brightening of the system visible as a hump 
around phase 1.0. Both systems may serve as ideal testers 
of mutual reflection treatment used in modern light curve 
synthesis programs. Indeed, the presence of small residuals 
visible near phase 1.1 in SC5 129441 and near the primary 
eclipse in SC6 67221 may (Fig. [^J suggest that simple re- 



flection used in the analysis (see §3.1) might be too crude a 
simplification. Unfortunately the use of the detailed reflec- 
tion (which should be more appropriate for these binaries) 
together with the eccentric orbits makes WD program prac- 
tically unsuitable for analysis due to an enormously large 
computational time (Wilson 1992). 

Both systems show no substantial third light contri- 
bution. The expected velocity semiamplitudes should be 
K 2 ^ 140 km/s for SC5 129441 and Ki ~ 150 km/s 
and K 2 ~ 140 km/s for SC6 67221. SC5 129441 was sug- 
gested by WW1 to have complete eclipses, but the present 
light curve analysis did not support this thesis. SC6 67221 
was not included by WW1. 

5.1.3 SC3 139376 

This binary seems to be the most massive and the most 
luminous in the sample (beside SC9 163575), consisting of 
two B0 V-IV stars of similar temperature: T ~ 32000 K. 
It is the only "classical" well-detached system included in 
the subset of likely distance indicators - the fractional radii 
of both components are below 0.2. The residuals of solu- 
tion are very small and only some larger scatter can be seen 
during eclipses due to very fast apsidal motion. The rate 
of the motion is ui ~ 5?8/year, so the apsidal period U is 
only ~ 60 years. Thus the analysis of this system can give 
parallel information about the components and an indepen- 
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dent check of stellar evolution models. The expected veloc- 
ity semiamplitudes should be K 1 ~ 190 km/s and K2 ~ 180 
km/s. 

5.14 SC6 215965 

The brightest binary in the sample (V — 13.94) with the 
spectrum BO V-IV + Bl IV. Table|2]suggests that the cooler 
component is more massive one and it is slightly evolved. 
The light curve solution showed that the contribution of the 
third light to this system might be neglected and photomet- 
ric parameters such as the temperature ratio (surface bright- 
ness ratio) could be very accurately determined. Also very 
small residuals suggested that the system had no essential 
intrinsic variability. This system was yet used by Harries 
et al. (2003) for a distance determination to SMC giving 
m — M = 18.83 ± 0.15. The physical parameters of both 
components are reasonably close to that presented in Ta- 
bic [J] It can be considered as an independent check of the 
method given in Section 2] and an additional argument for 
use this star for an accurate individual distance determi- 
nation. WW1 included this binary in their catalog on the 
basis of having complete eclipses but the present analysis 
(and results of Harries et al.) exludes this suggestion. 

5.2 SC5 38089 and the candidates suggested 
recently by Wyithe & Wilson 

5.2.1 SC5 38089 

Harries et al. (2003) determined absolute dimensions and a 
distance to this B0V+B0V system: m - M = 18.92 ± 0.19. 
Their light curve solution is essentially the same as the so- 
lution reported in this paper. However, there is an apparent 
discrepancy between the observational data on SC5 38089 
and the mass- luminosity relationship shown on Fig|H| Thus 
I would not recommended this star for individual distance 
determination, although it is very interested target for spec- 
troscopic investigations in the future. Table Ogives the esti- 
mate of physical parameters assuming that both components 
are normal main sequence stars. The expected B luminos- 
ity ratio is very close to unity what indeed was reported by 
Harries et al. 

5.2.2 SC4 103706 

Formally this binary also seems to have complete eclipses 
(it was included by WW1 in their subset of binaries showing 
totality) but the very large third light contribution (£3 ~ 0.4) 
may cause serious problems during modelling of the spectral 
energy distribution and thus I decided not to include this 
star into the most suitable candidates. 



5.2.3 SC4 163552 

WW2 included this binary according to the presence of com- 
plete eclipses, but as in the case of SC6 215965, the light 
curve analysis contradicted this thesis. The third light con- 
tribution to this system is sunstantial, at about 26%. This 
binary has small residuals, a very smooth light curve and 
can be considered as a "reserve" candidate. The mass ratio 



and temperature ratio are very close to unity. The expected 
velocity amplitude is K ~ 250 km/s for both components. 

5.2.4 SC6 11143 

WW1 classified this system as a well-detached. Indeed, my 
analysis showed that the components had the smallest frac- 
tional radii of all the sample (beside SC6 67221) and are 
normal B2 V stars. However, this binary happens to blend 
with another eclipsing binary (see §2.2) which contributes 
about 30% to the total light of system. Thus any detailed 
analysis of SC6 11143 should at the same time take into ac- 
count a careful analysis of the blending binary. Probably this 
includes simultaneous solution of multicolour light curves of 
both binaries what is, in principle, possible, but may give 
ambiguous results. 

5.2.5 SC6 319966 

This binary was included by WW1 according to the presence 
of complete eclipses which seems to be incorrect in the view 
of my analysis. This binary had the highest £3 contribution 
from all the sample - nearly 50 % - which is most probably 
too large for a detailed modelling of spectrophotometry. 



6 CONCLUDING REMARKS 

6.1 Third light contribution problem 

The observed difference of the £3 contribution to the to- 
tal light between binaries from A and B subsets may be 
understood as follows. The components of binaries from B 
subset are relatively much closer to each other than those 
from A subset and the proximity effects, depending on the 
gravitational distortion and overall geometry of the system, 
are much more distinct. Thus proximity effects give us ad- 
ditional constraints to the solution. When we consider the 
third light £3 we see that its influence on the light curve of 
well-detached systems is of little importance: equally good 
solutions can be obtained for wide range of £3 values. 

The situation is different for closer systems - small 
changes of £3 produce noticeably systematic O-C residuals, 
especially near quadratures. This is an important reason for 
preferring such systems over classical well-detached systems 
in determination of the distance using individual eclipsing 
binaries, especially in the case of the extragalactic bina- 
ries observed usually on a rich stellar background of a host 
galaxy. The probability of blending with another star(s) in 
such environment is high (crowding effect). Also the blend 
may be a member of the system itself as at least 30% of bi- 
naries (Batten 1973) are found to form hierarchical multiple 
systems of stars. Indeed, neglecting of the third light contri- 
bution in the case of HV 2274 and HV 982 may be one of the 
prime sources of small differences on the level 1.5<r between 
distances derived for three LMC systems as was announced 
by Ribas et al. (2002). 

6.2 Fast information about absolute dimensions 

The idea of deriving the absolute properties of binary stars 
from knowledge of their photometric properties alone has 
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long tradition (e.g. Gaposchkin 1938, Kopal 1959 - Chap- 
ter 7). The methods based on this idea use somewhat statis- 
tical approach: one or both components of a binary conform 
an empirical mass-luminosity relation and the temperature 
is imposed by observed spectral type. There are two impor- 
tant advantages of such statistical method 1) it may be used 
for large samples of binaries and 2) it gives fast, though ap- 
proximate physical parameters of binaries. Afterward we can 
choose from the sample most interesting systems and inves- 
tigate them more accurately using spectroscopy, spectropho- 
tometry or multicolor photometry. However, past methods 
neglected the interstellar reddening (which influences both 
the brightness and the colours) and, what is more important, 
the distance to the binary. As a matter of fact the distance is 
irrelevant to these methods and thus it can be derived sim- 
ply by comparison with the apparent magnitude m, even 
without any radial velocity measurements (e.g. Gaposchkin 
1968, Dworak 1974). Moreover, they based their reasoning 
on quite crude spectral type - temperature calibrations and, 
of course, as the temperature come in high powers into equa- 
tions involved in the method, the obtained parameters were 
usually very inaccurate. 

In present paper this "old-fashioned" method was in- 
novated to account for the interstellar reddening E(B — V) 
and recent B — V colour - temperature calibration (Flower 
1996). It turns out that such approach causes a big improve- 
ment in the temperature determination and overall physical 
parameters as was demonstrated by an accurate "reproduc- 
ing" of absolute dimensions of three LMC eclipsing bina- 
ries. However, it was done for a little bit larger distance to 
LMC than the distance derived, on average, from these LMC 
stars. A shift of the distance modulus A(m — M) = +0.1 
is quite small and comes probably from some inaccuracies 
of Flower's temperature or/and bolometric corrections cali- 
bration for hot stars. 

Anyway we should expect that this kind of method may 
give fast approximate determination of absolute dimensions 
for many binaries in the Magellanic Clouds, especially as we 
have large archives of binaries discovered by OGLE, MA- 
CHO and EROS projects at our disposal, and may be used 
also for eclipsing binaries in other close galaxies in the near 
future. 
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